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ABSTRACT 

The formation, merging, and accretion history of massive black holes along the hier- 
archical build-up of cosmic structures leaves a unique imprint on the background of 
gravitational waves at mHz frequencies. We study here, by means of dedicated simu- 
lations of black hole build-up, the possibility of constraining different models of black 
hole cosmic evolution using future gravitational wave space-borne missions, such as 
LISA. We consider two main scenarios for black hole formation, namely, one where 
seeds are light (~ 10 2 M©, remnant of Population III stars), and one where seeds are 
heavy (>10 4 M©, direct collapse). In all the models we have investigated, massive black 
hole binary coalescences do not produce a stochastic GW background, but rather, a 
set of individual resolved events. Detection of several hundreds merging events in a 3 
year LISA mission will be the sign of a heavy seed scenario with efficient formation 
of black hole seeds in a large fraction of high redshift halos. At the other extreme, a 
low event rate, about a few tens in 3 years, is peculiar of scenarios where cither the 
seeds are light, and many coalescences do not fall into the LISA band, or seeds are 
massive, but rare. In this case a decisive diagnostic is provided by the shape of the 
mass distribution of detected events. Light binaries (m < 10 4 M©) are predicted in a 
fairly large number in Population III remnant models, but are totally absent in direct 
collapse models. Finally, a further, helpful diagnostic of black hole formation models 
lies in the distribution of the mass ratios in binary coalescences. While heavy seed 
models predict that most of the detected events involve equal mass binaries, in the 
case of light seeds, mass ratios are equally distributed in the range 0.1 — 1 . 

Key words: black hole physics - cosmology: theory - early universe - gravitational 
waves 



1 INTRODUCTION 

Massive black hole (MBH) binaries (MBHBs) are among the 
primary candidate sources of gravitational waves (GWs) at 
mHz frequencies (see, e.g., Haehnelt 1994; Jaffe & Backer 
2003; Wyithe & Loeb 2003, Sesana et al. 2004, Sesana et al. 
2005), the range probed by the space-based Laser Interfer- 
ometer Space Antenna (LISA, Bender et al. 1994). Today, 
MBHs are ubiquitous in the nuclei of nearby galaxies (see, 
e.g., Magorrian et al. 1998). If MBHs were also common 
in the past (as implied by the notion that many distant 
galaxies harbour active nuclei for a short period of their 
life), and if their host galaxies experience multiple merg- 
ers during their lifetime, as dictated by popular cold dark 
matter (CDM) hierarchical cosmologies, then MBHBs in- 
evitably formed in large numbers during cosmic history. MB- 
HBs that are able to coalesce in less than a Hubble time (as 
defined at the epoch of their formation) give origin to the 



loudest GW signals in the Universe. Provided MBHBs do 
not "stall", their GW driven inspiral will then follow the 
merger of galaxies and protogalactic structures at high red- 
shifts. A low-frequency detector like LISA will be sensitive 
to GWs from coalescing binaries with total masses in the 
range 10 3 - 10 6 M out to z ~ 5 - 10 (Hughes 2002). Two 
outstanding questions are then how far up in the dark halo 
merger hierarchy do MBHs form, and whether stellar and/or 
gas dynamical processes can efficiently drive wide MBHBs 
to the GW emission stage. 

Today we know that MBHs must have been formed 
early in the history of the Universe. Indeed, the luminous 
z ~ 6 quasars discovered in the Sloan Digital Sky Sur- 
vey (Fan et al. 2001) imply that black holes more mas- 
sive than a few billion solar masses were already assembled 
when the universe was less than a billion years old. Sev- 
eral scenarios have been proposed for the seed MBH forma- 
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tion: seeds of m scc d ^fewxlOOM© can form as remnants of 
metal free (PopIII) stars at redshift J>20 (Volonteri, Haardt 
& Madau 2003, hereinafter VHM), while intermediate-mass 
seeds (m see( j ~ 10 Mg) can be the endproduct of the dy- 
namical instabilities arising in massive gaseous protogalac- 
tic disks in the redshift range 10<Jz<J15 (Koushiappas, Bul- 
lock & Dekel 2004, hereinafter KBD; Begelman, Volonteri 
& Rees 2006, hereinafter BVR, Lodato & Natarajan 2006). 
All these models have proved successful in reproducing the 
AGN optical luminosity function in a large redshift range 
(1 < ^ < 6), but result in different coalescence rates of MB- 
HBs, and hence in different GW backgrounds. 

In this paper we use the computational tools developed 
in Sesana et al. 2005, to characterize the expected GW sig- 
nal from inspiraling MBHBs in the different seed formation 
scenarios. Our aim is to understand the LISA capability to 
place constraints on MBH formation scenarios prior to the 
reionization epoch, looking for reliable diagnostics to dis- 
criminate between the different models. 

The paper is organized as follows. In § 2 we describe the 
different proposed seed formation scenarios. In § 3 we sum- 
marize the basic of the detection of GW from MBHBs. In 
§ 4 we compare the LISA detection rate and the properties 
of the detected MBHB population arising from the different 
seed formation scenarios. Finally, we summarize our main 
results in § 5. Unless otherwise stated, all results shown be- 
low refer to the currently favored ACDM world model with 
tt M = 0.3, D.A = 0.7, h = 0.7, fib = 0.045, a 8 = 0.93, and 
n = l. 



2 MODELS OF BLACK HOLE FORMATION 

In our hierarchical framework MBHs grow starting from pre- 
galactic seed MBHs formed at early times. The merger pro- 
cess would inevitably form a large number of MBHBs during 
cosmic history. Nuclear activity is triggered by halo merg- 
ers: in each major merger the more massive hole accretes 
gas until its mass scales with the fifth power of the circu- 
lar velocity of the host halo, normalized to reproduce the 
observed local correlation between MBH mass and veloc- 
ity dispersion (ttibh — cr* relation). Gas accretion onto the 
MBHs is assumed to occur at a fraction of the Eddington 
rate, as empirically shown in simulations of AGN feedback 
in a merger driven scenario. The scaling we adopt is based 
on the fitting formula by Hopkins et al. 2005 (see Volonteri, 
Salvaterra & Haardt 2006 for details). 

In this scenario there is a certain freedom in the choice 
of the seed masses, in the accretion prescription, and in the 
MBHB coalescence efficiency. 

In the VHM model, seed MBHs form with masses 
m scc d ~ fewxlO 2 M0, in halos collapsing at z = 20 from 
rare 3.5-cr peaks of the primordial density field (Madau & 
Rees 2001), and are thought to be the end-product of the 
first generation of stars. 

A different class of models assumes that MBH seeds 
form already massive. In the KBD model, seed MBHs form 
from the low angular momentum tail of material in halos 
with efficient molecular hydrogen gas cooling. MBHs with 
mass 
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We have fixed the free parameters in Eq. [T] by requiring 
an acceptable match with the luminosity function (LF) of 
quasars at z < 6. We note that, by requiring that the model 
reproduces the LF, the number of MBH seeds is very much 
reduced with respect to Koushiappas & Zentner (2006), 
where most of the black hole growth was due to black hole 
mergers. 

Here A is the so called spin parameter, which is a 
measure of the angular momentum of a dark matter halo 
A = J\E\ 1/2 /GM% 2 , where J, E and M h are the total an- 
gular momentum, energy and mass of the halo. The angu- 
lar momentum of galaxies is believed to have been acquired 
by tidal torques due to interactions with neighboring halos. 
The distribution of spin parameters found in numerical sim- 
ulations is well fit by a lognormal distribution in A sp in, with 
mean A sp m = 0.04 and standard deviation a\ = 0.5 (Bullock 
et al. 2001, van den Bosch et al. 2002). We have assumed 
that the MBH formation process proceeds until z m 15. 

In the BVR model, black hole seeds form in halos sub- 
ject to runaway gravitational instabilities. Gravitational in- 
stabilities are likely the most effective process for removing 
angular momentum. BVR have suggested that gas-rich ha- 
los with efficient cooling and low angular momentum (i.e., 
low spin parameter) are prone to global dynamical insta- 
bilities, the so-called "bars within bars" mechanism (Shlos- 
man, Frank & Begelman 1989). In metal-free halos with 
virial temperatures T v i r <;10 4 K, hydrogen atomic line emis- 
sion can cool the gas down to ~ 8000 K. In smaller halos, 
provided that molecular hydrogen cooling is efficient, gas 
can cool well below the virial temperature. The amount of 
material participating in the "bars within bars" instability, 
however, is much smaller in such mini-halos, leading to the 
accumulation in the proto-galaxy centre of only a few tens of 
solar masses. We assumed here, as in BVR, that MBH seed 
formation is efficient only in metal free halos with virial tem- 
peratures T v i r ^10 4 K. The "bars within bars" process pro- 
duces in the center of the halo a "quasistar" (QSS) with 
a very low specific entropy. When the QSS core collapses, 
it leads to a seed black hole of a few tens of solar masses. 
Accretion from the QSS envelope surrounding the collapsed 
core can however build up a substantial black hole mass very 
rapidly until it reaches a mass of the order of the "quasis- 
tar" itself, Mqss ^ 10 4 - 10 5 M Q . The black hole accretion 
rate adjusts so that the feedback energy flux equals the Ed- 
dington limit for the quasistar mass; thus, the black hole 
grows at a super-Eddington rate as long as Mqss > Mbh 
Mbh(<) ~ 4 x 10 5 (i/10 7 yr) 2 M i.e., M BH oc t 2 . 

In metal rich halos star formation becomes efficient, and 
depletes the gas inflow before the conditions for QSS (and 
MBH) formation are reached. BVR envisage that the pro- 
cess of MBH formation stops when gas is sufficiently metal 
enriched. Given the uncertainties in the efficiency of spread- 
ing metals, we consider here two scenarios, one in which 
star formation exerts a high level of feedback and ensures a 
rapid metal enrichment (BVRhf) and one in which feedback 
is milder and halos remain metal free for longer (BVRlf ) . In 
the former case MBH formation ceases at z m 18, in the lat- 
ter at z « 15. The BVRhf model appears to produce barely 
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enough MBHs to reproduce the observational constraints 
(ubiquity of MBHs in the local Universe, luminosity func- 
tion of quasars) . We consider it a very strong lower limit to 
the number of seeds that need to be formed in order to fit 
the observational constraints. 

Figure [1] shows the number of MBH binary coalescences 
per unit redshift per unit observed year, dN/dzdt, predicted 
by the five models we tested. Each panel shows the rates 
for different men = Tax + "12 mass intervals. The total coa- 
lescence rate spans almost two orders of magnitude ranging 
from ~ 3 yr^ 1 (BVRhf) to ~ 250 yr" 1 (KBD). As a general 
trend, coalescences of more massive MBHBs peak at lower 
redshifts (for all the models the coalescence peak in the case 
ibbh > 10 6 Mq is at z ~ 2). Note that there are no merging 
MBHBs with m B H < 10 4 M in the KBD and BVR models. 



2.1 MBHB dynamics 

During a galactic merger, the central MBHs initially share 
their fate with the host galaxy. The merging is driven by dy- 
namical friction, which has been shown to efficiently merge 
the galaxies and drive the MBHs in the central regions of 
the newly formed galaxy when the mass ratio of the satellite 
halo to the main halo is sufficiently large, that is when (to- 
tal) mass ratio of the progenitor halos is, P — M a /M>QA 
(Kazantzidis et al. 2005). 

The efficiency of dynamical friction decays when the 
MBHs get close and form a binary. In gas-poor systems, 
the subsequent evolution of the binary may be largely de- 
termined by three-body interactions with background stars 
(Bcgclman, Blandford & Rees 1980), leading to a long coa- 
lescence timescale. In gas rich high redshift halos, the orbital 
evolution of the central SMBH is likely dominated by dy- 
namical friction against the surrounding gaseous medium. 
The available simulations (Escala et al. 2004; Dotti et al. 
2006; Mayer et al. 2006) show that the binary can shrink to 
about parsec or slightly subparsec scale by dynamical fric- 
tion against the gas, depending on the gas thermodynamics. 
We have assumed here that, if a hard MBH binary is sur- 
rounded by an accretion disc, it coalesces instantaneously 
owing to interaction with the gas disc. If instead there is no 
gas readily available, the binary will be losing orbital energy 
to the stars, using the scheme described in Volonteri, Madau 
& Haardt (2003) and in Volonteri & Rees (2006). 



3 GRAVITATIONAL WAVE SIGNALS 

Full discussion of the GW signal produced by an inspiraling 
MBHB can be found in Sesana et al. 2005, along with all 
the relevant references. Here we just summarise the basic 
equations. 



3.1 Characteristic strain 

Consider a MBHB at (comoving) distance r{z). The strain 
amplitude (sky and polarization averaged) at the rest-frame 
frequency f r is (e.g., Thorne 1987) 
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Figure 1. Number of MBHB coalescences per observed year 
at 2 = 0, per unit redshift, in different »tibh = TO + m2 
mass intervals. Solid lines: VHM model; short-long dashed lines: 
KBD model; short-dashed lines: BVRlf model; long-dashed lines: 
BVRhf model. 



where M = m^m^/fmi + ni2) 1 ^ 5 is the "chirp mass" 
of the binary and all the other symbols have their standard 
meaning. The characteristic strain is defined as 
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where *Jn is the number of cycles spent in a frequency in- 
terval A/ ~ /. Equation 2] is valid only if the typical source 
shifting time (the time spent in a given A/ ~ / bin) is 
shorter than the instrumental operation time. This is al- 
most always the CELS6, cLS LISA would be sensitive to the 
signal emitted in the last 1 — 3 years before the MBHB co- 
alescence, and the operation time is expected to be >3 yrs. 



3.2 Resolved events 

An inspiraling binary is then detected if the signal-to-noise 
ratio (S/N) integrated over the observation is larger than 
the assumed threshold for detection. The integrated S/N is 
given by (e.g., Flanagan & Hughes 1998) 
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Here, / = /,-/(l + z) is the (observed) frequency emit- 
ted at time t = of the observation, and A/ is 
the (observed) frequency shift during the observational 
time t. Finally, h ruls is the effective rms noise of the 
instrument. The total LISA h Iuls noise is obtained by 
adding in quadrature the instrumental rms noise (given 
by e.g. the Larson's online sensitivity curve genera- 
tor |http://www.srLc altech.edu/~shane/sensitivity) and the 
confusion noise from unresolved galactic (Nelemans et al. 
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4.1 Event number counts 
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Figure 2. Rcdshift distribution of MBHBs resolved with S/N > 
5 by LISA in a 3-year mission. Line style as in Fig.[T] The number 
of events predicted by KBD model (long-short dashed curve) is 
divided by a factor of 10. The top-left corner label lists the total 
number of expected detections. 



2001) and extragalactic (Farmer & Phinney 2003) WD-WD 
binaries. Given the uncertainties on the very-low frequency 
LISA sensitivity, we adopt a pessimistic cut at 10 -4 Hz. We 
will discuss later the impact of changing the low-frequency 
cut-off of the sensitivity curve. 

Given a coalescence rate R, and the source frequency shift 
rate /, we can derive the number of individual binaries re- 
solved with S/N > s, i.e., (Sesana et al. 2005): 
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In equation [6] /isco is the observed frequency emitted at 
the Keplerian innermost stable circular orbit (ISCO), and 
/min <C /isco is the minimum observed frequency of the 
spiral-in phase. 



4 MBH FORMATION MODELS AND THE GW 
SIGNAL 

In this section we discuss the characteristics of the GW sig- 
nal detectable by LISA as predicted by the different models 
of MBH formation and evolution discussed in section 2. All 
the results shown here assume a LISA operation time of 3 
years, a cut-off at 10 -4 Hz in the instrumental sensitivity 
and a detection integrated threshold of S/N = 5 (eq. [5}. 



Figure [2] shows the redshift distribution of LISA MBHB de- 
tections. There are substantial differences between the differ- 
ent models. The KBD model results in a number of events 
(~ 700) that is more than an order of magnitude higher 
than that predicted by other models, with a skewed distribu- 
tion peaked at sensibly high redshift, z<;10. It is interesting 
to compare the number of detections with the total number 
of binary coalescences predicted by the different formation 
models. The KBD model produces ~ 750 coalescences, the 
VHM model ~ 250, and the two BVR models just few tens. 
A difference of a factor ~ 3 between the KBD and the VHM 
models in the total number of coalescences, results in a dif- 
ference of a factor of ~ 10 in the LISA detections, due to the 
different mass of the seed black holes. Almost all the KBD 
coalescences involve massive binaries (mi<;10 4 Mg), which 
are observable by LISA. The KBD and BVR models differ 
for the sheer number of MBHs. The halo mass threshold 
in the KBD model is well below (about 3 orders of mag- 
nitude) the BVR one, the latter requiring halos with virial 
temperature above 10 4 K. In a broader context, results per- 
taining to the KBD model describe the behaviour of families 
of models where efficient MBH formation can happen also 
in mini-halos where the source of cooling is molecular hy- 
drogen. 

It is difficult, on the basis of the redshift distributions 
of detected binaries only, to discriminate between heavy and 
light MBH seed scenarios. Although the VHM and BVRlf 
models predict a different number of observable sources, 
the uncertainties in the models are so high, that a differ- 
ence of a factor of two (96 for the VHM model, 44 for 
the BVRlf model) cannot be considered a safe discriminant. 
Moreover the redshift distributions are quite similar, peaked 
at z ~ 6 — 7 and without any particular feature in the shape. 
Arguably, the initial mass distribution of seeds is the main 
variable influencing the number of LISA events. However, a 
factor of two difference in the number counts, in fact, can be 
ascribed to uncertainties in other aspects of the MBH evo- 
lution, both dynamical and related to the accretion history. 
Sesana et al. (2005) showed how uncertainties in the dynam- 
ical evolution timescales reflects on the detectable events. 
Also the accretion history plays a role which may not be 
marginal. For instance, if MBHs grow fast (accreting around 
the Eddington rate or higher), the number of coalescences 
that can be detected by LISA is higher than, say, a case in 
which the accretion rate is sub-Eddington for most sources. 



4.2 Black hole masses and mass ratio distributions 

In Sesana et al. 2005 we showed that LISA will be sensi- 
tive to binaries with masses <J10 3 Mq up to redshift ten. 
Hence the discrimination between heavy and light MBH 
seed scenarios should be easy on the basis of the mass func- 
tion of detected binaries. This is shown in figure [3] As ex- 
pected, in the VHM model, the mass distribution extends to 
masses <J10 3 M©, giving a clear and unambiguous signature 
of a light seed scenario. VHM predict that many detections 
(about 50%) involve low mass binaries (m BH < 1O 4 M ) 
at high redshift (z > 8). These sources are observable dur- 
ing the inspiral phase, but their /isco is too high for LISA 
detection (see Sesana et al. 2005, figure 2). Heavy seed see- 
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Figure 3. Mass function of the more massive member of MBHBs 
resolved with S/N > 5 by LISA in a 3-year mission. Line style 
as in figure [T] All curves are normalized such as the integral in 
d log(mi ) gives the number of detected events. 
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Figure 4. Mass ratio distribution of MBHBs. long-dashed curve: 
< z < 3; short-dashed curve: 3 < z < 10; dotted curve: z > 10; 
so/id curve: all redshifts. Le/f panels: VHM model, all coalescences 
(upper panel), and coalescences detectable by LISA with S/N > 5 
in a 3-year mission (lower panel). Most high-redshift events with 
mass ratios of order unity involve light binaries which cannot 
be observed by LISA. Right panel: KBD model; almost all the 
coalescences can be observed with a S/N > 5. 



narios predict instead that the GW emission at /isco, and 
the subsequent plunge are always observable for all binaries. 

In figure [4] we show the mass ratio distribution of the 
resolved events for the VHM and the KBD models. Model 
KBD (as well as BVRlf and BVRhf, not shown here) pre- 
dict a monotonically increasing distribution, the majority of 
detections having a mass ratio q — M%jM\ ^ 1. A large 
fraction of observable coalescences, in fact, involve MBHBs 
at z > 10, when MBHs had no time to accrete much mass 
yet. As most seeds form with similar mass (~ 1O 4 M0, see 
KBD; ~ 10 2 M , see VHM), mergers at early times involve 
MBHBs with q ~ 1. In massive seeds scenarios, almost all 
coalescences are observable, and the mass ratio distribution 
is dominated by z > 10 mergers between seeds (q ~ 1). In 
scenarios based on Population III remnants, z > 10 merg- 
ers involve MBHs with mass below the LISA threshold. The 
detectable events happen at later times, when MBHs have 
already experienced a great deal of mass growth. VHM mod- 
els therefore produce a mass ratio distribution which is flat 
or features a broad peak at q ~ 0.1 — 0.2, depending on the 
details of the accretion prescription. This is due to the fact 
that both the probability of halo mergers (because of the 
steep DM halo mass function) and the dynamical friction 
timescale increase with decreasing halo mass ratio. Hence, 
equal mass mergers that lead to efficient binary formation 
within short timescales (i.e., shorter than the Hubble time) 
are rare, while in more common unequal mass mergers it 
takes longer than an Hubble time to drag the satellite hole 
to the center. 

All the results shown above assume 3 yrs observation 
and a cut-off in the LISA sensitivity curve at 10 -4 Hz. Even 
in the pessimistic case of a 1 yr mission lifetime, however, 
with a sensitivity curve cut at 10 -4 Hz and assuming the 
BVRhf seed model, LISA is expected to observe at least 
two or three MBHB merging events. We stress once again 
that the BVRhf model provides a strong lower limit to the 
number and redshift distribution of forming seeds, on the 
basis of current observational constraints. 



4.3 Confusion noise 

If the number of merging sources is so large that there are, on 
average, at least eight sources above threshold per frequency 
resolution bin, then the total signal will be observed as a con- 
fusion noise (Cornish 2003). A detectable confusion noise of 
cosmic origin would provide much information on the emit- 
ting population, but, on the other hand, would be added (in 
quadrature) to the instrumental noise, reducing the interfer- 
ometer capability of detecting individual sources. Assuming 
a mission lifetime of three yrs, the predicted confusion noise 
(see Sesana et al. 2005 for details), varies by an order of 
magnitude for different models, but lies, for all models, be- 
low the LISA sensitivity curve. In the pessimistic view of a 
one year mission, the confusion noise is enhanced roughly by 
a factor of three. As shown in figure [5] the confusion noise 
predicted by the KBD model is expected to be comparable 
to the rms noise at frequencies <C3 x 10~~ 4 Hz. If the sen- 
sitivity curve cuts-off at 10 -5 Hz; the quadrature addition 
of such a noise would result in a slight decrease of the total 
LISA sensitivity in the frequency range 3 x 10~ 5 - 3 x 10~ 4 
Hz. 
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Figure 5. Predicted confusion noises assuming a LISA op- 
crating time of one year. Thin lines: LISA rms confu- 
sion noise (solid line), as the quadratic sum of the LISA 
instrumental single-arm Michelson noise (dotted line, from 
http://www.srl.caltech.edu/~shane/sensitivity), and the confu- 
sion noise from unresolved galactic (Nelemans et al. 2001, long- 
dashed line), and extragalactic (Farmer & Phinney 2003, short- 
dashed line) WD-WD binaries. Thick lines: predicted confusion 
noises for the different MBHB models wc tested. Line style as in 
in figure [T1 



5 SUMMARY AND CONCLUSIONS 

Using dedicated Montecarlo simulations of the hierarchical 
assembly of DM halos along the cosmic history, we have 
computed the expected gravitational wave signal from the 
evolving population of massive black hole binaries. The im- 
print of black hole mergers and coalescences on the LISA 
data stream depends on the specific assumptions regarding 
MBH formation, and on the recipes employed for the hole 
mass growth via merger and gas accretion. 

We have considered two main frameworks for MBH for- 
mation, namely, one where seeds are light (~ 1O 2 M ), and 
one where seeds are heavy (£10 4 M o ). In the former, MBH 
seeds form at z ~ 20 with masses of few hundreds solar 
masses, and are though to be the endpoint of the evolution 
of metal-free massive stars (VHM). In the heavy seeds sce- 
narios, MBHs form in the centers of high-redshift gas-rich 
halos where angular momentum losses are efficient. KBD ex- 
plore a model where angular momentum is shed via turbu- 
lent viscosity, in all halos with efficient molecular hydrogen 
cooling. This seed formation scenario is very efficient, and 
predicts that seeds are widespread, forming in halos as small 
as a few 10 s Mq , provided that the total angular momen- 
tum of the halo is small enough. BVR explore a different 
scenario where angular momentum is transported via run- 
away gravitational instabilities ("bars-within-bars"). BVR 
envisage that the process would be more effective in halos 
with efficient atomic cooling, that is with virial temperature 
Tvh->10 4 K, and mass M h >10 8 Mq. MBH seeds are therefore 



much rarer in the BVR model. BVR envisage that the pro- 
cess of MBH formation stops when gas is sufficiently metal 
enriched. Given the uncertainties in the efficiency in spread- 
ing metals, we consider here two scenarios, one in which star 
formation exerts a high level of feedback and ensures a rapid 
metal enrichment (BVRhf), one in which feedback is milder 
and halos remain metal free for longer (BVRlf). In the for- 
mer case MBH formation ceases at z ~ 18, in the latter at 
z « 15. 

We have shown that, in all considered models, MBHB 
coalescences do not produce a stochastic GW background, 
but rather, a set of individual resolved events. A large frac- 
tion (depending on models) of coalescences will be directly 
observable by LISA, and on the basis of the detection rate, 
constraints can be put on the MBH formation process. De- 
tection of several hundreds events in 3 years will be the sign 
of efficient formation of heavy MBH seeds in a large fraction 
of high redshift halos (KBD). 

At the other extreme, a low event rate, about few tens 
in 3 years, is peculiar in scenarios where either the seeds 
are light, and many coalescences do not fall into the LISA 
band, or seeds are massive, but rare, as envisioned by, e.g., 
BVR (see also Lodato & Natarajan). In this case a decisive 
diagnostic is provided by the mass distribution of detected 
events. In the light seed scenario, the mass distribution of 
observed binaries extend to ~ 10 3 Mq , while there are no 
sources with mass below 10 4 Mq in the high seed scenario. 
Finally, we have shown that a further, helpful diagnostic 
of MBH models lies in the distribution of the mass ratios 
in binary coalescences. While heavy seed models predicted 
that most of the detected events involve equal mass binaries, 
in the case of light seeds, mass ratios are equally distributed 
in the mass ratio range 0.1 — 1. 

Should the early black hole population be dominated by 
massive systems (e.g., KBD), the GW signal can be accom- 
panied by an electromagnetic counterpart (EM, Miloslavl- 
jevic & Phinney 2005, Dotti et al. 2006, Kocsis et al. 
2006), in principle detectable by future high-sensitivity X- 
ray telescopes (e.g., XEU^I). The identification of an EM 
counterpart would have crucial implications for cosmol- 
ogy. LISA sources could then be used as "standard sirens" 
(Schutz 1986) to estimate fundamental cosmological param- 
eters (Schutz 2002; Holz & Hughes 2005), providing con- 
straints complementary to cosmic microwave background 
experiments, and supernovae search experiments. 

We note, however, that from the astrophysical point 
of view, even in absence of EM counterparts, by adopting 
the standard Wilkinson Microwave Anisotropy Probe cos- 
mological parameters (Spergel et al. 2003) the redshift of 
the merging MBHs will be known with the same precision 
as most sources observed in electromagnetic bands. 

In conclusion, concerning the detection of low frequency 
gravitational waves, massive black hole binaries are certainly 
one of the major target for a mission as LISA. On the astro- 
physical ground, LISA will be a unique probe of the forma- 
tion, accretion and merger of MBHs along the entire cosmic 
history of galactic structures. 

We would like to thank A. Vecchio for helpful discus- 
1 www.rssd.esa.int/index.php?project=XEUS 
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